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Fermions constitute an important component of matter and their quantization in presence of dy- 
namical gravity is essential for any theory of quantum gravity. We revisit the classical formulation 
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form for gravity together with the Nieh-Yan topological term which keeps the nature of Barbero- 
Immirzi parameter independent of inclusion of arbitrary matter with arbitrary couplings. With 
dynamical gravity, a priori, there are two distinct notions of 'parity'— orientation reversing dif- 
feomorphisms and improper Lorentz rotations. The invariance properties of the action and the 
canonical framework are different with respect to these and gravitational origin of parity violation 
seems ambiguous. 
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I. INTRODUCTION 



Incorporation of fermions in the background free quantum theory gravity has been dis- 
cussed in the hterature [ij. In the formulation of general relativity in terms of real SU(2) 
connection, Thiemann discussed loop quantization of standard model fields The fermions 
were treated in the second order form i.e. fermions couple to gravity through the spin con- 
nection (torsion free Lorentz connection). Perez and Rovelli returned to fermions in presence 
of the Hoist term and found that its coefficient playing the role of the inverse of the Barbero- 

n n 

Immirzi parameter, 7, becomes classically observable |3[. Mercuri discovered that with a 
further addition of suitable non-minimal fermionic couplings, 7 can be made classically un- 
observable. He also noted that these added terms (Hoist plus non-minimal) can be expressed 
as the Nieh-Yan topological term once the connection equations of motion are used. The 
strategy of adding non-minimal couplings to keep 7 classically unobservable was followed for 
N = 1, 2 and 4 supergravities also [5]. Canonical analysis and loop quantization of fermions 
with non-minimal couplings was discussed by Bojowald and Das [6]. It was subsequently 
realised that 7 will automatically be classically unobservable provided it is the (inverse) of 
the coefficient of the Nieh-Yan term (a total divergence) in the Lagrangian density. Thus, 
instead of the Hoist terms alone, if the Nieh-Yan term (Hoist + torsion^ piece) is used in 
conjunction with the Hilbert-Palatini, then for arbitrary matter and their couplings, the 7 
will drop out of the classical equations of motion. Furthermore, it is possible to systemati- 
cally derive the real SU(2) Hamiltonian formulation from such an action Q]. Fermions were 
also included in the canonical analysis. A necessary condition for a topological origin of 7 
is thus satisfied. The canonical analysis leading to real SU(2) formulation has since been 
extended to supergravities jsl] as well as further inclusion of the other two topological terms 
namely the Pontryagin and the Euler classes [o], [l^ . 

Fermions are also tied with possible parity violations jo], There are two distinct no- 
tions of 'parity': one related to orientation of the space-time manifold and one related to 
the improper Lorentz transformation. Depending upon the definitions of the basic canon- 
ical variables (with or without the sgn(e) factors in this work (section |V])), the canonical 
framework and the action are (non-) invariant under one of the notions of parity. These pos- 
sibilities are not distinguished in the previous works. The work in [3] is a little incomplete 
in the constraint analysis although the final results are correct. The constraint expressions 
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are also not in a form which is suitable for loop quantization. This work seeks to fill these 
gaps. 

We re-derive the real SU(2) formulation including a Dirac fermion. The analysis is done 
using restricted fields corresponding to 'time gauge'. Elimination of the second class con- 
straints, leads to the usual formulation. When a Dirac fermion is included, the solution of 
the second class constraint leads to non-trivial Dirac brackets between the SU(2) connection 
and the fermions. One can however make natural shifts in the definition of the connection to 
recover the canonical brackets. This also simplifies the constraints. Four Fermi interaction 
terms however survive in the Hamiltonian and are signatures of first order formulation. In 
the second order formulation where fermions couple to the torsion free connection, there are 
no terms quartic in the fermions. 

The straightforward derivation also introduces factors of sgn(e) := sign(det(e^)) (which 
equals iVsign(det(V^*)) in the time gauge parametrization) in appropriate places. For in- 
stance, we are naturally lead to the definitions: := sgn{e) y/qV^" , := sgn(e) +77 
(eqn. ( IT4l) ). Under an improper orthogonal transformation (0(3)) acting on the index i, the 
triad changes sign and so does the sgn(e) factor leaving invariant. This is as it should be 
since the index i on E represents adjoint representation of SO (3) while on V it represents 
the defining representation. For S0(3), both are equivalent but not for 0(3). Under an 
inversion, A*^ = — quantities in the defining representation change sign while those in 
the adjoint don't. The same reasoning applies to the definition of the connection. Now the 
connection is also even under inversion. The sgn(e) factor also change the behaviour of E^ 
and under the action of orientation reversing diffeomorphisms. These factors of sgn(e) 
however occur only in the intermediate derivations, the final form of constraints and basic 
variables are independent of these factors. 

The paper is organized as follows. Section [TTl deals with the Hamiltonian formulation for 
purely gravitational sector starting with the Hilbert-Palatini plus the Nieh-Yan action. Here 
steps in the analysis are given and a gap in the constraint analysis in [t] is filled in. The 
constraints are presented in the more standard form. The sign factors which appear in the 
action itself, are followed through in all expressions. Section IIIII deals with inclusion of a 
single Dirac fermion. Constraint analysis as well as simplification of constraints is given here. 
The section ends with the final form of the canonical formulation in the standard notation. 
Section |V] discusses the two distinct notions of parities and the invariance properties of the 
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action and the canonical framework. Section IVT] contains brief concluding remarks. 



II. THE HAMILTONIAN FORMULATION 

The starting point is a choice of (tensor /spinor) fields and a corresponding generally 
covariant, local action on 4 dimensional space-time M ~ [R x S3. The next step is to carry 
out a 3 + 1 decomposition to identify the Lagrangian which is a function of (tensor) fields 
on E3 together with their velocities with respect to the chosen time coordinate. The fields 
whose velocities appear in the Lagrangian are potentially the configuration space variables 
while those without velocities appearing in the Lagrangian are Lagrange multipliers whose 
coefficients will be primary constraints. This Lagrangian leads to the kinematical phase 
space. Now a constraint analysis a la Dirac is performed. If there are second class constraints, 
one may hope to simplify the analysis by solving the second class constraints. However, now 
one must use the Dirac brackets. These may not have the canonical form for the remaining 
variables (i.e. may not be Darboux coordinates) and a new choice of variables may be 
necessary. This is particularly relevant for Lagrangians which are linear in the velocities 
such as the Hilbert-Palatini-Nieh-Yan and the Dirac Lagrangians which typically do have 
primary, second class constraints. The classical Hamiltonian formulation is completed when 
the action is expressed in the Hamiltonian form together with first class constraints. We also 
have fields coordinatizing the kinematical phase space (after the second class constraints are 
eliminated) with the configuration space coordinates identified. 

A. Pure gravity 

We begin with the Lagrangian 4-forms built from the basic fields the co-tetrad e^^dx^ and 
the Lorentz connection ui^'dx^. 



Cupie,uj) = — sgn{e)-£ijKLR^-\co) Ae'^ Ae^ 



(1) 




(2) 
(3) 
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Here R :— diu + 00 Aou and T :— de + ou Ae are the usual curvature and torsion 2-forms^ and 
K, — SttG. The factor of sgn(e) is present because only then the Hilbert- Palatini Lagrangian 
matches with the y^\R{g). This arises from noting that determinant of the co-tetrad is 

given by, e =sgn(e)y^[g|. 

A 3 + 1 decomposition can be expressed using, 

d := d± + d\\ := dtdt + dx"'da 

= ei + e{ , u^-^ = (^/■^ + cuy^'^ etc (4) 

which leads to, 

>Cgrav = ^ md A d± {sgn{e)u}\\ij - r]ujpj) + 2r77|J| A d^e\\i + 2776^ A £>||7>|||| + 
uJ±i.j A [D\\ (sgn(e)E{f - T^Sjf) - 2i^e\ A 

+ A |e|[ A ^sgn(e)^7j|||| — | + surface terms] where, (5) 

E-'^"^ := A e"^, D := d + cuA is the Lorentz covariant derivative and e is the determinant 
of the co-tetrad, assumed to be non-zero. The ~ denotes the Lorentz dual defined for any 
quantity antisymmetric in a pair of Lorentz indices, X^"^ , as. 
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A 3-1-1 decomposition is carried out as usual by choosing a foliation defined by a time 
function T : M — )■ IR and a vector field t^(9^, transversal to its leaves. The vector field is 
normalised by t ■ dT = 1 so that the parameters of its integral curves, serve as the time 
coordinate. Given such a decomposition, we choose a parametrization of the tetrad and the 
co-tetrad as, 

el = Nn' + N^V^ , = , n'rij = -1 , nV/r^.j = ; (6) 

with V^V^ = Si + mn-^ , V;V,,' = 5^ . (7) 



^ Our conventions are such that (i) r] = 7"^, our action matches with Ashtekar-Lewandowski (AL) action 
including signs and factors. Our identification of K^^ is same is that of AL while our F* is minus that of 

AL. To match with the standard notation we first change A — > rjA and then set 77 = 7"^. This results in 
standard definition of our — > 7"^^/*^^, the diffeomorphism and the Hamiltonian constraints remain 
unchanged. The symplectic form becomes "j~^EfdtA\. Therefore the Poisson brackets also become the 
standard ones. 
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Expressing the wedge products in terms of the components, using dx* Adx^ Adx^ Adx'^ :— 
£tabc^^t ^^1 ^^2 g^^^ Separating the IJ sums into (Oi), (jk) sums, we write the 

gravitational Lagrangian density as, 

^tabc 

>Cgrav = ^5— [2V^V^dt (sgn(e)a)„oi - V^^aOi) + V^V^dt {sgn{e)uaij - V^aij) 
+UJto^ {Da (sgn(e)SO: - r^E^l) - (sgn(e)n° - V^t) } 

{d^ (sgn(e)Etl - v^}) } 
+N |n° (p;* (sgn(e)/2o^6c - vRoibc)) + 

(sgn(e)^oi6c - vRoibc^ + {sgn{e)Rijbc - -Roibc)) } 
{V^ (sgnie)Roibc - vRoibc)) + 
K' (sgn(e)^o*c - vRoib^j + (&gTi{e)Riji^ - r^-Ri^bc)) } 

The terms in the last two lines, come from the torsion piece of the Nieh-Yan term. We 
have also defined. 

We will now restrict to configurations such that Ui — 0, no = —1. This also implies that 
= = Vq and that are invertible with as the inverse^. We also define the 3- 
metric qab '■= Va^b^ij (which is positive definite in classical theory) and denote q :=dct{qab) ■ 
Many terms in the above equation drop out. In particular, there isn't any time derivative 
of Vajo in the last square bracket. We also promote the torsion components to new 
independent variables T"* and to V^, for the terms in the last square bracket. The number 
of variables is restored back by introducing two primary constraints: :—V^ — V^^O and 

:= f"- — T"' with ^f, 0^ as the corresponding Lagrange multiphers. Note that T"° is 
not promoted to a new variable. This simplifies the Lagrangian to, 

^tabc 



^ This would correspond to the choice of the so-caUed time gauge if we started without restricting the 
configurations a priori. 
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-rutoi {Da (sgn(e)E°: - r^E^:) - (sgn(e)n° - rjT^^ 
+u:u,{Da (sgn(e)ni-r;E;i)} 

|k (sgn(e)^oi6c - vRoib^ } 

+7V«' |v;;,\/i- (sgn(e)i?,,,e - r^i?,,,,) } 

K - 

+N { {DaT'^') } + 7V«' {v:, {Dam ]] (10) 

At this stage it is convenient to introduce new notations for certain combinations of the 
components of the Lorentz connection as well as those of the (co)tetrad. These are, 



Et e{elet-e\el) ^ = sgn{e)^qVt (11) 

K , l^'jk^a'' (12) 

'■^aij — ^ijk^a > '^aOi — ' '^aij = ^ijlK^ (13) 



and we have used e = A'"det(yj) = A'"sgn(e)y^. Note that N > implies that sgn(e) = 
sgn(det(\4*)) and therefore we continue to use sgn(e) every where. 

With these, the first velocity term becomes, 2E'^dt (sgn(e)ii'* — r/F^) suggesting the iden- 
tification of a canonical pair, (^4^, E!^) with 

Al := sgn{e)Kl-r^Tl (14) 

The velocity terms become, 

(>Cgrav)velocity = ^ dtA^ + d^V^ (15) 

At this stage we have 10 Lagrange multipher fields, N, N"", u/"^ and 18 canonical variables, 
Al^, Ef while the co-tetrad and the Lorentz connection constitute 13 -|- 24 = 37 fields. Nine 
more fields are yet 'unaccounted'. If we define 

uj' := sgn(e)a;/-^ + rjuj' , 

we can see that, = 
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The nine remaining variables, MkiX^y ^ire identified through the definition^, 



sgn e 



{ElC-EiC) + -S'^''MuEl , Mu = Mi, , E:E': = 6t , E^E^ = 6] . (16) 



The 18 components of uj"^ are organized as the 9 A]^, the 3 C ^i-nd the 6 Mm- The components 
K, r of the Lorentz connection u can be expressed in terms of the A, (andM using the inverse 
formula. 



ij 



P, I J cIJ KL 

KL 



1 ( ( \^ IJ V 

iT^ (^Sn(e)u;. - -< 
The K, r are given in terms of the independent canonical variables A, (, M as, 

Y^2{'Sn{e)Al + 'l{M^^Ei + sgn{e)S'r,Eie)} , = \sgn{e)K:-Al). (17) 

For later convenience, we eliminate F in favour of K, A. We get. 



grav 



1 



E^dtAi + TtdtK 



KV, 



a 



^ a ' a ^ 



Kg'' 

Khia' 

kU 



rjia rpa 

-'-i -'-i 1 



2 



Da ( sgn{e)tll - r/S°M - ( sgn(e)Sf, - r/S^° 



be 



- VIT 



irpaO 



kij 



s 



tabc 



a' 



£ 



tabc 



^tabc 

( sgn(e)i?oifec - vRoibc 



DaT 



aO 



(18) 
(19) 
(20) 

(21) 

(22) 
(23) 
(24) 



In anticipation, we will refer to the Q^'^ as the boost constraints, Qk as the rotation con- 
straints, Ha' as the diffeomorphism constraints and H as the Hamiltonian constraint although 
at this stage the interpretation of the transformations generated by these is ambiguous. The 
other two terms in the second line are the two additional primary constraints which identify 
the canonical coordinates T, V with the torsion and the co-tetrad components respectively. 



^ It is possible to bypass this decomposition and work directly with and as independent canonical 
variables. The conjugate momenta of Afy , will be subsumed by the conjugate momenta of /iT* (ic|. We 
will however continue to use these variables. 
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We have also used Um := Aj, uuj := £ijk^^ ■ Straightforward algebra then leads to^ 



Qk 



irpaO 



sgn{e)daE''' - sgn{e)C - V^T 



Ha' := r,F:,,El + 



7] 



X' {daft - Ka{rS + V''£^,kf^K^Z<e)Kl - Al)] where, 



Kb 



H 



(25) 
(26) 

(27) 
(28) 

(29) 



As there are no velocities of Mij^C^ variables, we have their conjugate momenta vr*-' 
~ TTfc as additional primary constraints. In Dirac's terminology, the total Hamiltonian will 
be the minus terms in the second line of eqn 0181) plus linear combinations of vTj ^ TCij ^ 0. 
All these are to be preserved during an evolution and this leads to secondary constraints 
and determination of some Lagrange multipliers. 

Consider the preservation of the constraints vTj, iiij, v^, rf ^ 0. This leads to the equations: 



{TTfc, Htot} ~ 
{7lkl,Htot} ~0 



A,- 



5C 
60' 



Oi 



+ A 



+ N 



5H 



(30) 
(31) 
(32) 
(33) 



Now we note that = sgn(e)5^ and {vl{x) , Tj (y)} = 5* 5^ 5^{x,y) (1 + r]'^)^^. This 
means that equation ( 130|1 can be solved for Aj, equation (132|) can be solved for 0^ and 
equation ( !33|) can be solved for (This also means that {ni,Q^^) and (t;*,rj') constraints 
form second class pairs.) Furthermore, from the appendix, we note that the boost and the 
rotation constraints are independent of M'^' and therefore the first two terms of eqn ( l3T|l 
are zero. This equation will turn out to give a secondary constraint. These statements 



We set K = 1 for notational convenience. It will be restored back later. 
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remain true also in the presence of fermionic matter. Let us use the abbreviation H :- 
J N'^'Ha + Nl-i. Then, using equation fl32|) . we obtain, 

l + r/2 

Now, 



j A.<J« + XgX - H] (34) 



It is easy to see that for X = j K-Q + \- Q, the Poisson bracket with t;* is weakly zero, with 
or without matter (fermionic). Using these, the equation ( 13T1) becomes, 

sgn(e) ^ _ ^^.^y.'^) and (36) 



^M'^' 2(1 + r/2) 

^ = 9.(iV''K) + £:Vn^Ar^r^ + NKl so that, (37) 

H (38) 



5 r/ sgn(e) . 5 
2— + ^'"'^^^ ^ 

7] sgn(e' 



The second hne of the last equation is independent of matter contribution. These appear 
only in the first line of that equation. Noting that the M^' dependence in H appears only 
through K*, we can trade derivative w.r.t. M*'' to that with respect to K^^. Furthermore, 
the A^, dependence is both explicit as well as implicit through K]^. This simphfies equation 
to, 

r fx 

{E,,5\ + E,,5i) J- ^ {5,,V,- + 5uV^ - 25,iVt) {da{N'Vi) + S^,,ViN'V^^ + NK\} (39) 

On the left hand side, the 5 signifies that only the explicit dependence on A\ is to be picked 
up. The right hand side is independent of matter contributions. The sgn factors cancel out. 
A somewhat lengthy but straight forward computation yields, 



^ ^^ + — {EakQi + EaiGk) where, (40) 

Ski 1 



. ^ {sgn{e){Mki - 5kiM\) - v{AakE^ + A^iE'^ - 26kiAi^E1)} 
rj 1 + rj"^ 

+ £'""'{VakdbVa + Vaid,Vck) (41) 
10 



All the dependence on the derivatives of the lapse and shift variables disappears. The 
secondary constraint is just Ski ^ 0: 

= TT^ {sgn(e)M^^- - {A^,E^^ + A^E^) - 6'^ (sgn(e)M^, - 2r^A^Ef ) } 

+ vS""'^ {V%Vl + Vid,V^:) (42) 
sgn(e)M*^- = -{l + r]')S'''''{v:d,V^ + V^d,V:-6'^V,%V,k}+v{AlE'^' +AiE'^') (43) 

When fermionic matter is added, this constraint gets modified and is discussed in the next 
section. 

There are no further tertiary constraints. 

From eqn. (!30l) . we notice that (vrfc,^°') forms a second class pair of constraints, one of 
which is a canonical variable. Defining a Dirac bracket relative to this pair will allow us to 
set these constraints strongly equal to zero. Because vTfc is a canonical variable, the Dirac 
brackets of variables other than vr^, C', among themselves, coincide with the corresponding 
Poisson brackets. This follows by noting (schematically), 

A := {n,G{q\p„n,C)} 
W,P„vr,C), 9{q\Pj,7rX)h ■= {f,9} + {f,n}A-'{G,g}-{f,G}A-'{n,g} 

^ U\Pj}* = {q\Pi} , (44) 

Since there is no explicit vr^ dependence in any of the remaining constraints, these ex- 
pressions remain the same. We impose = strongly and eliminate C = daE"-^ — V^*T°". 
From eqn. (IMTD of the appendix, this sets T''' = I^^^Y/. 

Similarly, {7Tki,S^^) form a second class pair hence we can define Dirac brackets relative 
to these and impose these strongly. The constraints rcki being canonical variables, the Dirac 
brackets among the remaining canonical variables remain the same as their Poisson brackets. 
Once again, setting TCij = changes no expressions but S^^ = gives M*-' in terms of Al_ and 
Ef. Also S'^ = implies that T'" = and sets the = f^ ^ 0. 

But now we can define Dirac brackets relative to the (i;* , rj*) second class pair. In this r° 
is now a canonical variable and therefore the Dirac brackets among the remaining variables, 
A^, E'^, remain equal to their Poisson brackets. 

When fermionic matter is included, the first two steps of elimination of second class 
constraints will remain the same but with non-zero torsion and this will lead to non-trivial 
Dirac brackets among the final set of variables in the third step. 
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Note: Solving Sij = and = strongly, determines Mij and in terms of A, E. 
Substitution in the expression for (eqn. [17]), leads to, 

cbcd pi 

O „ , jkj-^jr. j-,bk 



Using the identities, S'^'^E^ = S'^'^'V^'Vj and ^^''(l^^^V'/ - V'^^I^^^) = ^^''r^^f *;„„V;™V;" , one 
can see that this is precisely the the usual expression for in terms of the triad alone, 

K = {d,V^, - + V^d^V,,] . (45) 

Thus, in the pure gravity case, we get the Hamiltonian formulation in terms of the 
A\, Ej variables with the original Poisson brackets, we are left with only the rotation, the 
diffeomorphism and the Hamiltonian constraints and the total Hamiltonian is made up of 
these alone. It remains to simplify the expressions for these constraints. 

B. Simplification of Constraints 

After the second class constraints are imposed, we get T^°' = — sgn(e) rjVf'Q^ ~ , 
Tf = ^ and V^* = V^. The boost constraint, eqn. (j25l) . is strongly zero and the 
rotation constraint, eqn ( l26l) . takes the usual form: 

g' = T]daE''' + £'^^AiE''\ (46) 

Consider the diffeomorphism constraint. The second term in equation fl27|) simplifies as, 

( l±jr\ .ri sgn(e) Kg^El , 

= \^y^ — 7^ — (-^ ^^•^^') 



2 



[r^e + ^'^'A^^E^) 



1 + r/^ 



1 + 11^ 

V:,K,,V^ 6;^sgn(e) (47) 



r] 

.-. -Ha' = vK'bE'l - -V:,K,,V,'g' sgn(e) ^ vF^,,E^ (48) 
V 



In the last but one line we have used the equation ( ]A23p . 
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Next consider the Hamiltonian constraint (eqn.(l29|)). after using the expressions for tor- 
sion. 



n := -- 



1 E'^.E] 



3 k 



2 



£'iFi} - da (sgn(e) vVrQ') + 



l_E}El 
2 



-sgn(e) {vimi - d^Ki) + S'^MbK - A^K^)] 



(49) 



The /T-dependent terms can be written as (after taking out the (1 +rf )/{2rf) factor) and 



using 



= sgn(e)^''-K;, 



— 2sgn(e) 7]£-'\di 



EjE^i 



+ 2sgn(e) r]Ki£''''nVa 



al 



- 2 



sgn(ej 



{A\En£,,,8'^SKlE'^) + 



sgn(ej 



(KE'^nSr^kS'mniKlE 



The total derivative term simphfies as, 



[£^''A,,E^ + r]C^] = 



E]EIK[ 



sgn(e) 
1 + r/2 

96(sgn(e) V^g' 



sgn(e) 



(50) 



Using the equation for Mj^ and the boost constraint for ehminating C*, we get, 

sgn{e)KiE^^ - A'^E'^^ = -"Lg-^'^ {2V^^d,V^ - 5'^V^d,V,,] . 

Writing Aj^E"'" = {A^^E""" - sgn{e) Ki^E""") + sgn{e) K'^E""" in the 3rd term, the algebraic 
K— dependence simplifies to, 

-E'^E^ [KIK^, - KlKt) . 
Combining all the terms, the Hamiltonian constraint takes the simplified form. 



n := - 



1 E]El 



2 



7]- 



+ -da {sgn{e)Vrg') . (51) 



The total derivative term differs from eqn. 2.24 of 12|, because that expression is derived 
from the Hoist action while ours is derived from the Hilbert-Palatini-Nieh-Yan action. 

With this we recover the usual form of the constraints (in the time gauge) starting from 
the Hilbert-Palatini action with Nieh-Yan terms added. In the next sub-section we add a 
Dirac fermion minimally coupled to gravity. 
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III. ADDITION OF A DIRAC FERMION 



The Lagrangian we begin with is, 



r - --\ \ 

•'-Dirac 2 



Ae^7'^/.(c^, A, . . .)A - D^{uj, A, . . .)Ae^7'A (52) 
1 



D^{uj)\ := df,\ + -uJ-^aijX + ie'A^X + ...X , 
DiM^ ■■= I^^At + ^ujJ'X^crlj - te'A,,X^ + . . . A+j 7° 



(53) 



The . . . refer to possible couphngs of the Dirac fermion to other gauge fields eg the Maxwell 
field. These are suppressed in the following. The conventions for the 7^— matrices (space- 
time independent) are given in the appendix lA 1[ The factor in front of the Dirac Lagrangian 
is —i because our metric signature is (- + + +). This is crucial for the Dirac brackets and 
subsequent passage to quantization via 'Quantum brackets = ih Dirac brackets' rule, with 
the quantum brackets being realised on a Hilbert space. 

For future convenience we introduce \& := q^^'^X, \E't := g^^^A^. This absorbs away the 
factors in the Lagrangian as well as in the constraints. Note that the terms involving 
the derivatives of cancel out. The A fermionic variables being of density weight zero, the 

fermionic variables are of density weight 1/2. From now on we will use the half density 
variables. 

Substituting the 3 + 1 parametrization of the tetrad and using the time-gauge, the 
Lagrangian can be written as, 

/:Dirac = ^(^^5t^ - 5i(^t)^ _ ^uujG'/ - N^'u^, - NUp where, 

Q'^ = Q , Qf = - ^S^.k^a^"^ (54) 

= -^(^7°^a'^-:D:^7°^) , = '-V^i^YD^^-D^^Y"^) (55) 



As in the vacuum case, the fermionic contribution to the boost and the rotation con- 
straints is independent of K and hence of M and (. Hence, the Q^^ pair of second class 
constraint can be eliminated in exactly the same manner as before leading to the same de- 
termination of ( and of course without affecting the Poisson brackets among the remaining 
variables. Furthermore, as in the vacuum case, the Aj, A*, terms drop out from the equation 
( 13T|) . The secondary constraint is thus determined from equation ( 139|) with H now including 
the fermionic contribution. We only need to pick out the explicit A dependence. 
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This is easily done and leads to, 



^ J , (56) 

^ (i^..^f + - '-^^ {5.^h.^) (57) 



Adding these to the eqn. fl5^ gives, 



[5h - sgn(e) r/4^^t^5*] + — [^afc^^^* + ^a;6^f ] ~ (58) 



Hence the secondary constraint becomes. 



S,, ^ sgn(e) r75,,^t^^^ sf. (59) 
and the solution for Mjj becomes, 

sgn(e) M^^' = -(1 + 7]^)S'''" {V%Vi + Vld^V^ - b'^ld^V,^ + ^ {A\E''^ + A^.E'^^) 

_(ljtZZ!l sgn(e) 5*^'^t^5^ . (60) 

After imposing the second class constraints strongly, we get, 

ft = Tt = \S.,V^^ = sgn(e) v^^^l^"^ 

and this leads to non-trivial Dirac brackets. While we could simplify the remaining con- 
straints as before, it turns out to be better to check the Dirac brackets relative to the 
(i)* , ) second class pair which suggests a new definition of A and simplifies the constraints 
considerably. 

For the fermions, the action being linear in velocities, we have primary constraints, vr^ ~ 
A^, TTxt ~ A, which are second class. Also these variables fail to be Darboux coordinates - do 
not have vanishing Poisson brackets^ with the gravitational variables due to the ^/q factor. 
The shift to variables makes the matter and gravitational variables Poisson-commute. 

Defining Dirac brackets relative to these primary, second class constraints allows us to use 



strictly Generalised Poisson brackets [l3| . due to the Grassmann nature of the fermions. 
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as basic variables with Dirac brackets given by^, 

{vl/"(x), ^l{y)U = -^S''^S\x,y) (61) 
Now we define Dirac brackets relative to the {v, r) constraints: 

< = K - Vl rt = ft- ^S^.V"^ =ft- '-^^^Vr^i^^j,^ . 

The Dirac bracket is defined as, 

{f{x),giy)r := {f{x),g{y)}-jdz{f{x),vl{z)}{Tt{z),g[y)} 

+ j dz{f{x),Tt{z)}{vl{z),g{y)} . (62) 

Due to the presence of the fermionic term in rf , the Dirac bracket between A and the 
fermions is non-trivial while all other basic Dirac brackets remain the same as the corre- 
sponding Poisson brackets. Specifically, 

{Al(x),vl/"(y)}* = -sgn(e) l^El{^,^T5\x,y) , 

{A^,{x)MMY = +sgn(e) '-^E:{^^^,)j\x,y) . (63) 

This suggests that if we use, 

a; := Al-sgn(e) |i?:v[;t^5^ , 

all Dirac brackets become standard Poisson brackets. Thus, while simplifying the constraints 
we will first effect the shift A = A+ ^^^^^^ ^ E'^^'jr,'^. This results in very natural expressions 
for the constraints. 

The solution f l60|) for M*-' simplifies as, 

M'\A,E,<if) = M!^{A,E)-^6'^<if^-f5'^ where, M(A. E) is defined in (Ii3l) . (64) 
while the K, T simplify as, 

K{A,E,C) = KliAEX) , K{A,EX) = rM,EX)-'-^El^^^,^ (65) 

^ The minus sign in the basic Dirac bracket, is correlated with the sign in the Lagrangian. Upon quan- 
tization, we will have anti-commutator, [^',^'^]_|_ ~ +h which is consistent with the Hilbert space inner 
product. 
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The underlined functions are the functions obtained in the pure gravity case ('torsion- free'). 
In particular, £^(^4, E, () = TiV) as given in equation fHSj) . when the second class constraints 
are strongly solved. This is relevant for regularization of the Hamiltonian constraint later 
on. 

The shift of A does not affect the expression for the boost constraint (1A12I) since the 
extra term is killed by the £. For the same reason the rotation constraint (]A13|) is also 
unaffected. 

The diffeomorphism constraint simplifies as follows. Recall eqn.f l27p . 

V 

-^{^7°D„,^-D:;^7°^} (66) 

where, 

= sgn(e) V.'^vl/t^^vl/ , T'^o= -rj sgn{e)Vtgi^, ^ sgn{e)vVrg'p , G'^ = ^h^aoi^ 
Da'"^ = (^da'^ - '-A^^, 75^0.^) + K, + sgn(e)^^ ao^^ 

D:;^ = {da:^ + ^A^, ^^75^0.) - K, ^ (^1 + sgn(e)^^ ^o, (67) 

In the above expression, we have separated the /T— dependent terms and also solved the 
second class constraints. 

The first step in simplification is to put A\ = + sgn(e)^£'*\E'^75\E'. The K independent 
pieces just have the A dependence being replaced by A dependence with one additional 
term, 

[ir-indep. terms] (A) = [fsT-indep. terms] (A) - sgn(e)|E*,((9a£;f )(^"^75^) (68) 

The i^'— dependent terms combine to give, 

[X-dependent terms] = - sgn(e)^E*,^f (^^75^) (69) 

In getting this simplification, we have used the boost constraint, the rotation constraint as 
well as = K{A, M, () expressions. 

The middle term of the first line of equation f l66|) and the last term of equation fl68|) . com- 
bine to give — sgn(e)|i?*/\E'^75\E'^j^'^'^. This in turn combines with the simplified if— dependent 
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term to give a term proportional to the total rotation constraint which is weakly zero. Thus, 
finally, we get the simplified diffeomorphism constraint as, 

-Ha'iA, E, ^) = vFUA) E\ - {^7°^^a'^ - ^7°^} where, (70) 

Va'^ := {da>^ -i^l^ao^^^ , V:;^ := {d,,^ - i ^ ^^^a^^j (71) 

Notice that this is a simple additive form of the vacuum diffeomorphism and the Dirac 
diffeomorphism constraint (coupled to the SU(2) gauge connection A). Four fermion terms 
generated by the second class constraint have been neatly absorbed in the shifted A. 
Lastly, we simplify the Hamiltonian constraint. Recall eqn.f l29|) . 



-sgn{e)£^\ [^{d,K[ - d^Ki) + S'^MTK - A^^K^)) } + (72) 
+ '-Vt {^I'D^m - :D:^y^} (73) 

As before, second class constraints have been solved and the i^'— dependent terms are sepa- 
rated. Substitute A\ = ^ + sgn(e)2E*^t^^^. Thanks to eqn. (165|) . we replace everywhere 
K ^ I£. The A— dependent terms generate additional terms. 

2 ^ ^ ihA^) ^ 



2 



'-^'\fUA) + ^-^V%E^i + ^1 sgn(e) + ^ {¥^,^y 



7 H 



and. 
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The dependent terns in the fermionic Hamiltonian and the torsion pieces combine to 
give, 

-sgn(e) S^^K^E'^^^ + sgn(e) {vVtO^F) + ^^^a^^" (^^75*) 
The terms generated by the shift in A, simphfy to, 

2r/^ 2r]^ ^ 16 ^ 

Using Qp — Q^^ — Q^^^ and dropping the term containing , leads to 



1 E'^El 



2 Vq 



c f 1 -I- -r)2 



16 ^ 



(74) 



+ 



1 4* /^'* 



^H5^)sgn(e) + 



Q'pQkF 



AVq " 2^ r/ 

The first two hnes of the above equation are exactly the same as the Hamiltonian for 
vacuum case and therefore are simplified in exactly the same way as before. In particular, 
the derivatives of X terms after partial integration, generates another total derivative term, 
namely, (r/~"^ + ?7)(9a(sgn(e) V^"^*^^), which combines with the last term in the first line leaving 
us with ri~^da{sgn{e)V^Ql^^) . The /c— dependent terms just produce minus the second term 
in the first line as before. 

The third line is the contribution from the Dirac Hamiltonian coupled to the A field while 
the terms in the fourth line are the extra term including 4-fermions terms. 
Using the relation. 



4^ " 2^ V 

The final form of the Hamiltonian constraint is. 



n 



1 E'^^E' 



2 ^/g 



l + rf 



{KiKH - KiKt) + -da (sgn(e) ¥^01,,) 



V 



-KIE', 



3 gkg^^ 



16 ^ 



2r?^/g 

, where g'p = ^^^75^0^* 



(75) 
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The Hamiltonian constraint thus consists of additive combination of the vacuum and 
the Dirac Hamiltonian (coupled to A). However, unlike the diffeomorphism constraint (rfO|) . 
there are the additional terms in the square bracket. These extra terms contain contributions 
quartic in the fermions as well as quadratic in fermions. Notice that there are no explicit 
factors of sgn(e) in the final expressions. The appears which is the same as in the vacuum 
case and therefore the properties needed in using the Thiemann identities in the quantization 
of the gravitational Hamiltonian constraint continue to hold. We have also checked that if 
the connection equations of motion for the spatial components are solved and the solution 
for these torsion components are substituted back in the action, the above Hamiltonian is 
recovered. 

We close this section by reverting to the standard notation eg 12 1. 



For this, we first substitute — )■ rjAl^ (and — )■ K^). This removes the factors of t]~^ 
from the definition of given in eqn. ( 128|1 . It also extracts a common factor of f] from the 
rotation constraint, Q^. The diffeomorphism and the Hamiltonian constraints remains the 
same except for removing the factor of rj^^ from the total derivative term in the Hamiltonian 
constraint. The symplectic term, E'^dtA''^ — )■ {riE^)dtA''^. And finally we put rj = 7"^ and 
restore k, (Recall the overall k,~^ factor in the Lagrangian density of eqn f|T8l) ). 

Here are the final expressions: 



■[7i,7fc]; 



(76) 

(77) 
(78) 

(79) 



n := {£'\fL{A) - (1 + 7^) {KlKl - KiKl) } + 7^. (sgn(e) YtQ^.) 



KPA (M/t^^M; 



3 (^Ws*) 

-K- 
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(^t^^OrO^^)(^t^^^O.^) 



(80) 



Note: Dimensionally, k L\ {A,K,d) ^ L-\ E ^ , P ~ L-^ , (^, ^) 
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The inverse square root of q and appearing above are manipulated exactly as in the 
vacuum case. As remarked earlier, the A and K above correspond to the vacuum case for 
which the Thiemann identities hold. Explicitly, the identities we would use are: 

sgn(e)£:''-K = 5^ ' ? (det(K))' = det(£;f) (81) 



\/d<Ef) 



E^E'^ of /"I 

S'^'' , ' = —S'^^lAiix), / d^y^} . (82) 

HE{i) ---^j ^^'i^^c ^ ^ J ^S^i^)KPr (83) 
K = {k^')-' !^He{1) J AVg} , sgn(e)X:(x) = ^(a;) , (84) 

These identities suffice to derive a quantization the Hamiltonian constraint from that of the 
'Euclidean Hamiltonian constraint' (the first term in the Hamiltonian constraint) and of the 
volume operator. 

IV. ACTION OF CONSTRAINTS AND THEIR ALGEBRA 

It is easy to see that the gauge constraint generates correct gauge transformation of the 
basic fields. Specifically, with ^(A) := j^^d^xh^Qi, 

Ka(x),e;(A)} = -P„A^ = - daK^ - Sy,A^,K' (85) 

{P^{x),G{m = ^if^'n (86) 

{^"(x),^(A)} = -iA*(75ao.^)" (87) 

{¥^{x)M^)] = +^A*(*^5<7o^)a (88) 

If we compute the infinitesimal action of the Ha constraint on the basic variables, we see 
that it equals the Lie derivatives of the basic variables only up to an SU(2) gauge transfor- 
mation. We are however free to modify the constraints by adding suitable combinations of 
themselves. So we define the diffeomorphism constraint as: 

C{N) := / d^xNXa with, 

Ca := Ha - KQi = Pida^, " ^^(^l^') + \ (^^^a* " ^a*^ ' *) (89) 
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which leads to the infinitesimal transformations, 





),C(iV)} 


— /■ - 4* 

— /^]^J-^^ 


= daiN'A^) + N\dbA',-daAl) 


(90) 




),C(iV)} 


= C -P" 


= N^dbPt - PidbN" + 1 ■ {dbN^)P^ 


(91) 




),C(iV)} 


= ^A^^" 




(92) 




),C(iV)} 




= iV^9feVl/t + i . (9,iV^)vi>t 


(93) 



This implies that {var, J N°'Ca} = C]\]a{var) for all variables. The Gauge constraint 
already generates the correct gauge transformation of the basic variables. By inspection, 
it follows that the gauge constraints (weakly) commute with the diffeomorphism and the 
Hamiltonian constraint, the gauge constraint and the diffeomorphism constraints form sub- 
algebras and the diffeomorphism constraint transforms the Hamiltonian constraint by the 
Lie derivative. The non-trivial bracket is the bracket of two Hamiltonian constraints. 

V. PARITY AND INTERNAL PARITY 

Recall that we begin with the (co-)tetrad field e^, the Lorentz connection uj^^^ and the 
fermion fields A, A (or \E', ^) defined over a manifold M ~ IR x S3 which is assumed to 
be orientable. With the topology specified, M can be taken to be time-orientable with 
respect to all the metric tensors constructed by the parametrization (!6|7|) . Obviously, S3 is 
orientable as well. 

There are two distinct sets of discrete transformations: orientation reversing diffeomor- 
phism of M and a 0(1,3) transformation with determinant = -1. We will keep the time 
orientation fixed. Orientation reversing diffeomorphism of M will then be reversing the 
orientation of S3. We will refer to these as parity transformations. The improper Lorentz 
transformations A^, will also be taken so that detA = — 1 and A°q = 1 and will be referred 
to as an Lorentz parity transformation. 

After going to the canonical framework in the 'time gauge', we have the fields A^, P", ^ 
defined on S3. The parity transformations are the orientation reversing diffeomorphism 
of S3 and the improper 0(3) transformation, inversion (say) will be the Lorentz parity 
transformation. 

In the Lagrangian framework, the Hilbert-Palatini action (also the Euler and the cos- 
mological terms) are invariant under both sets of transformations while the Nieh-Yan (as 

22 



well as the Pontryagin and the Hoist) actions change signs under parity but are invariant 
under Lorentz parity. Hence the combined action is invariant under Lorentz parity and 
non-invariant under parity. 

The variables of the canonical framework are defined in terms of those of the Lagrangian 
framework. These definitions of the SU{2) connection in terms of K and F and the conju- 
gate momentum in terms of the triad are consistent with the 5*0 (3) gauge transformation 
extended to include the Lorentz parity. Thus the triad which transforms by the defining 
representation changes sign under Lorentz parity. The 'densitized triad' (or the conjugate 
momentum) transforms by the adjoint representation and should be invariant under Lorentz 
parity. The sgn(e) factor in their definitions precisely takes care of this. The same can be 
seen in the definition of the connection. It is easy to see that the symplectic structure and 
the constraints (vacuum) are all invariant under Lorentz parity. 

When fermions are included, these are scalars under diffeomorphism and transform as 
^ — > 7°^ under Lorentz parity. All the constraints including fermions are invariant under 
Lorentz parity. This is true in both the Lagrangian and the canonical frameworks. 

With regards to parity the situation is different. The action is not invariant under par- 
ity, due to the Nich-Yan term. In the canonical framework, the connection is not simply 
even/odd under parity since the K term changes sign while the F does not. The 'densitized 
triad' also acquires an extra minus sign under parity (behaves as a 'pseudo- vector of weight 
1'). The symplectic structure thus is not invariant. The constraints also are not invariant 
under parity. This is consistent with the non-invariance of the action. 

The action is invariant under parity combined with 7 —7. Our definitions have the 
appropriate factors of 7 to restore the simple (even) behaviour of the basic canonical variable 
resulting also in the invariance of the Poisson brackets and constraints. 

In short, Lorentz parity is an invariance of the action as well as the canonical framework 
and parity is not. However, parity combined with 7 — > —7 is an invariance of both action 
and the canonical framework. It is not our intention to suggest that the combined operation 
be a physical symmetry (which depends on the quantum theory) , but it is useful in checking 
the algebra. 

One could try to change the definitions of the basic variables by dropping the sgn(e) 
factors. This will result in expressions which can be obtained from the above by putting 
sgn(e) = 1. This will restore the 'densitized triad' to its usual density weight 1 vector density 
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status and the connection to its 1-form status. The canonical framework is then invariant 
under parity (without changing sign of 7). The action however is still not invariant. 

Under Lorentz parity, the connection does not have simple behaviour and the conjugate 
momentum will be odd. This results in non-invariance of the canonical framework. The 
action however is still invariant under Lorentz parity. If the sign of 7 is changed along 
with the Lorentz parity transformation, then the basic variables are even, the symplectic 
structure is invariant and so are the constraints and the action. Thus, the definitions without 
the sgn(e) factors, interchanges the role of Lorentz parity and parity appropriately combined 
with 7 — )■ —7. 

With two different definitions (with and without the sgn(e) factors), we can ensure either 
invariance of the canonical framework under parity or under Lorentz parity and this is 
independent of minimally coupled fermions. The action however unambiguously remains 
invariant under Lorentz parity and non-invariant under parity. Which of these is more 
appropriate? 

Observe that if we were to consider formulation in terms of the metric tensor, then the 
notion of Lorentz parity is not even definable as there is no internal Lorentz transformation. 
On the other hand, existence of fermions (spinorial fields) requires an orientable manifold and 
using the tetrad formulation making the Lorentz parity notion available. If the orientation 
of the manifold is regarded as a fixed background structure, then parity transformation is 
excluded by definition and Lorentz parity alone is available. Which of these is relevant from 
an observational point of view is not very clear and so also the issue of 'parity violation' via 
gravitational interactions. 



VI. CONCLUDING REMARKS 



The next natural step is the loop quantization of the fermions with interacting gravity. 
The kinematic Hilbert space of this system has been already constructed |2;]. The procedure 
given by Thiemann can be followed in toto. The extra feature, not available in Thiemann's 
discussion are the quartic terms in the fermionic sector. Their regularization has been given 
in jo] and we don't have any thing new to add to this. 

In summary, we have presented a canonical form of a Dirac fermion, minimally coupled to 
the tetrad form of gravity including the Nieh-Yan term. The Canonical analysis shows that 
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the coefficient of the Nieh-Yan term is the inverse of the usual Barbero-Immirzi parameter 
while the topological nature of the Nieh-Yan term guarantees its non-appearance in the 
classical equations of motion. One could consider additional non-minimal couplings, not 
affecting the symplectic structure, but these will not change the status of the Barbero- 
Immirzi parameter. We saw the natural appearance of the sgn(e) factors which also serve 
to give definitions consistent with the two distinct notions of parity and Lorentz parity. 
There is no Lorentz parity violation either in the Lagrangian framework or the Hamiltonian 
framework for this system. As noted earlier, among matter fields, fermions alone are sensitive 
to the SU(2) action and contribute to the rotation (Gauss) constraint. Its implications for 
the homogeneous, diagonal models has already been explored by Bojowald and Das [6]. For 
the same reasons, fermions are likely mediators in the black hole evaporation process. Due to 
the role of fermions in the chiral anomalies in the usual quantization, they are also probes to 
see how loop quantization does or does not accommodate chiral anomalies. For exploration 
and elaboration of these issues, it is necessary to have a sufficiently precise control over the 
fermion-gravity system and classical analysis is the first step in this direction. 
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Appendix A 



1. Notation and Conventions 

The basic fields are the co-tetrad e^, the Lorentz connection a;^^ and a Dirac spinor A°. 
The co-tetrad is taken to be non-degenerate and defines a metric gf^u :— ej^e'^rju, rju — 
diag(— 1, 1, 1, 1). The general coordinate indices, ji^v . . are raised/lowered with g^j,i,, g'^'^ 



while the Lorentz indices, I, J, ■ ■ ., are raised/lowered with rju, r]^'^ . 

e det(ei) := ^S^'^-^SuKLeiete^e^^ , S'^'^ = 1 = £:oi23 ; (Al) 

2eJ:Z e{e^^e^j-e^^e^j) = ^S^^'^^SuKLe^e^^ (A3) 

X" := ^^"^iX^^ , X = -X (A4) 

D^ei - D^^ , D,ei - d^e'^ + uj^'jci , (A5) 

277^1 = 7V + 7V , \ ' 75 := ^7^S'7' (A6) 



A := At70 , 7) = 707,7° , = 2^75^0, , (A7) 

D»A := a^A + ^a;/VjA , 'DiM^ := j^^At + la;/^AVl^+| (A8) 

An explicit representation for the Dirac matrices is chosen to be. 




2. Torsion Components 

The torsion is defined geometrically as the covariant derivative of the co-tetrad. These 
involve the components of the Lorentz connection which have been obtained as functions of 
the canonical variables A, E, M, (. Similarly the boost and the rotation constraint are also 
expressed as functions of A,E,(- Eliminating ( and some combinations of S,A,E, we can 
express the torsion in terms of the constraints and the V, T variables. Subsequently, we will 
set y = y and T = T to simphfy the expressions for the torsion as well as the constraints. 
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The equations are, 



rj-iai 

r" 



be 



V cabcrpi 

2 



Oi 
vac 



-•vac 



Substitution of, 



^abc^k 



^aiEj 



= sgn{e){daE^'-C)-V:T-'' 
= TidaEl-SujKE'^' +£hjKT''' 

sgn(e)A,iE; + |Mi,- + %gn{e)£ijkC' 



l + r/2 
1 



-r|Aa^E'; + 



sgn(ej 



(AlO) 
(All) 
(A12) 
(A13) 

(A14) 
(A15) 

(A16) 



leads to, (using E°; = sgn(e)ygl^j" and i?* = sgn(e)\4Vyg) 



aO 



T 



l + r/2 



V;'^sgn(e) {r^C + A,,El] 



-r^Vt (sgn(e)^:,, - r^^c) + V sgn{e)£^ ,,VtV^t 



bk 



(A17) 



l + rf 



{sgn(e) {M'^ - 5'm\) - {A'^E'^^ + A^.E'^^ - 25^^A^E,") } 



'a a 



d^E^'' + 



7] 



ck /\m rpbn 
2^ mn^b ^ 



[See eqn.dUD] 



(A18) 



1 + r7 

7] 



7]^ 



[sgn{e)E\V^T'^' + £% [r^Ql^ + sgn(e)Cc)] 



£1 \ v£'mnT"'V- + sgn(e)K'T'^° + vGL + sgn(e)^?; 



^vac 



rpaiyj 



\S^^ + |sgn(e)^^iCc 



(A19) 
(A20) 



In the last hne of the equation! ] Al 91) we have used equation ( 1A17I) . Now if we further put 
T = T,V = V, T""^ simphfies to give, 

T^' = - rj sgn{e)VrgL, , T^^ = ^S'^v; + sgn{e)^£%V;g'^',, (A21) 

Using these in the expressions of the boost (1A12P and the Gauss constraints (lAlSp . lead to. 



->0i 

'vac 



sgn{e) 



daE''' - 



c 

1 + 



+ 



V 



V 



_pi A3 p 



<ak 



(A22) 
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W = Y^,C + Y^^S'.kA^aE^' (A23) 
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